Symmetry and Supersymmetry in Nuclear Physics 

A. B. Balantekin 

Department of Physics, University of Wisconsin, Madison, WI 53706, USA 



Summary. — A survey of algebraic approaches to various problems in nuclear 
physics is given. Examples are chosen from pairing of many-nucleon systems, nu- 
clear structure, fusion reactions below the Coulomb barrier, and supernova neutrino 
physics to illustrate the utility of group-theoretical and related algebraic methods 
in nuclear physics. 

PACS 03 . 65 . Fd - Algebraic methods. 

PACS 21.60.Fw - Models based on group theory. 

PACS 25 . 70 . J j - Fusion and fusion-fission reactions. 

PACS 26 . 30 . -k - Nucleosynthesis in novae, supernovae, and other explosive envi- 
ronments. 



1. Introduction 

Symmetry concepts play a very important role in all of physics. Originally symmetries 
utilized in physics did not change the particle statistics: such symmetries either transform 
bosons into bosons or fermions into fermions. Natural mathematical tools to explore 
symmetries are Lie algebras, i.e. the set of operators closing under commutation relations, 
schematically shown as 

(1) [Gb, Gb] = Gb- 
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Lie algebras are associated with Lie groups. One could think of Lie groups as the expo- 
nentiation of Lie algebras. 

On the other hand, supersymmctrics transform bosons into bosons, fermions into 
fermions, and bosons into fermions and vice versa. Natural mathematical tools to explore 
them are superalgebras and supergroups: Superalgebras, sets containing bosonic (Gb) 
as well as fermionic (Gf) operators, close under commutation and anticommutation 
relations as shown below: 

[Gb,Gb] = Gb, 
[Gb,Gf] = Gf, 

(2) {Gf,Gf} = G b 

The simplest superalgebra can be easily worked out. Consider three dimensional 
harmonic oscillator creation and annihilation operators and define 



(3) 



K o = l(b\bi + l) K + = \bU = (K^. 



It is easy to show that the operators defined in Eq. (3) satisfy 

(4) [K , K±] = ±K±, [K+, K_] = -2K . 

The algebra depicted in Eq. (4) is the SU(1,1) algebra [1]. Casimir operators are opera- 
tors that satisfy the condition 

(5) [C(G), K ] = = [C(G), K±] 

Casimir operators obtained by multiplying one, two, three elements of the algebra are 
called linear, quadratic, cubic Casimir operators. For SU(1,1) the quadratic Casimir op. 
is 

(6) C 2 = Kl-\ (K + K_ + K_K + ) 

The concept of dynamical symmetries found many applications in nuclear physics (see 
e.g. [2]). Consider a chain of algebras (or associated groups): 

(7) G 1 dG 2 d---dG 71 

If a given Hamiltonian can be written in terms of the Casimir operators of the algebras 
in this chain, then such a Hamiltonian is said to possess a dynamical symmetry: 



(8) 



n 

ff = ^[a i C 1 (G i ) + ftC 2 (G i )] 

i=l 
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Obviously all these Casimir operators commute with each other. Consequently, find- 
ing the energy eigenvalues of this Hamiltonian reduces to the problem of reading the 
eigenvalues of the Casimir operators off the existing tabulations. 

Introducing spin (i.e., fcrmionic) degrees of freedom represented by the Pauli matrices 
as well as the bosonic (harmonic oscillator) ones and defining 

(9) f + = i£y4, F_ = i£>A, 

i i 

one can write the following additional commutation and anticommutation relations: 

[K , F±] = ± l -F ± , 

[K+, F + ] = 0=[K_, F_], 

[K±, F T ] = T F±, 
{F±, F ± } = K ± , 
{F + , F_} = K 

The operators K + , K_, K , F + , and F_ generate the Osp(l/2) superalgebra, which 
is non-compact [1, 3]. Since the operators K + , K_, K alone generate the SP(2) ~ 
SU(1,1) subalgebra we can write the group chain 

(10) Osp(l/2) D SU(l, 1) D SO(2). 

(Note that the operator Kg can be viewed as the Casimir operator of the SO (2) subalgebra 
of SU(1,1)). The Casimir operators of Osp(l/2) and Sp(2) - SU(1,1) are given by 
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(11) c 2 (Osp(l/2)) = i(L+|) 2 = ij 

(12) C 2 (Sp(2)) = iL 2 - A, 

where L is the angular momentum carried by the oscillator, i.e. Lj = ti^jVk = i e ijkb\.bj. 
It can easily be shown that a harmonic oscillator Hamiltonian with a constant spin-orbit 
coupling 

(13) H= 1 -(p" + ^)+X^-L+ 3 - 
can be written in terms of the Casimir operators of the group chain given in Eq. (10): 

(14) H = 4\C 2 (Osp(l/2)) - 4AC 2 (Sp(2)) + 2K . 



This provides perhaps the simplest example of a dynamical supersymmetry. 
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2. Fermion pairing 

Pairing is a salient property of multi-fermion systems and as such it has a long his- 
tory in nuclear physics (for a recent review see [4]). Already in 1950 Meyer suggested 
that short-range attractive nucleon-nucleon interaction yields nuclear ground states with 
angular momentum zero [5]. Mean field calculations with effective interactions describe 
many nuclear properties, however they cannot provide a complete solution of the under- 
lying complex many-body problem. For example, after many years of investigations, we 
now that the structure of low-lying collective states in medium-heavy to heavy nuclei 
are determined by pairing correlations with L=0 and L=2. This was exploited by many 
successful models of nuclear structure such as the Interacting Boson Model of Arima and 
Iachello [6, 7, 8]. 

Pairing plays a significant role not only in finite nuclei, but also in nuclear matter and 
can directly effect related observables. For example, we know that neutron superfluidity 
is present in the crust and the inner part of a neutron star. Pairing could significantly 
effect the thermal evolution of the neutron star by suppressing neutrino (and possibly 
exotics such as axions) emission [9]. 

Charge symmetry implies that interactions between two protons and two neutrons 
are very similar; hence proton-proton and neutron-neutron pairs play a similar role in 
nuclei. In addition isospin symmetry implies that proton-neutron interaction is also very 
similar to the proton-proton and neutron-neutron interactions. Currently there is very 
little experimental information about neutron-proton pairing in heavier nuclei as one 
needs to study proton rich nuclei to achieve this goal. But one expects that data from 
the current and future radiative beam facilities will change this picture. 

First microscopic theory of pairing was the Bardeen, Cooper, Schriffcr (BCS) theory 
[10]. Soon after its introduction, the BCS theory was applied to nuclear structure [11, 
12, 13]. Application of the BCS theory to nuclear structure has a main drawback: BCS 
wave function is not an eigenstate of the number operator. Several solutions were offered 
to remedy this shortcoming such as adding Random-Phase Approximation (RPA) to the 
BCS theory [14], projection of the particle number after variation [15], or projection of 
the particle number before the variation. The last technique was recently much utilized 
in nuclei (see e.g. [16]). Pairing correlations may also play an interesting role in halo 
nuclei [17]. It should also be noted that the theory of pairing in nuclear physics has many 
parallels with the theory of ultrasmall metallic grains in condensed matter physics (see 
e.g. [18]). 

2T. Quasi-Spin Algebra. - The concept of seniority was introduced by Racah to aid the 
classification of atomic spectra [19]. Seniority quantum number is basically the number of 
unpaired particles in the j n configuration. Seniority-conserving pairing interactions are 
a very limited class, however such interactions make an interesting case study. Kerman 
introduced the quasi-spin scheme to treat such cases [20]. In this scheme nucleons are 
placed at time-reversed states \j m) and (— l)^~ m ^\j — m). Introducing the creation and 
annihilation operators for nucleons at level j, a, m and a,j m , the quasi-spin operators are 
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written as 



(15) S+=^(-l)^a} m a]_ m , 

m>0 



(16) Sr =Y,(-l) {j - m) a j - m a jm , 

m>0 

and 

( 17 ) Sj = \ { a J rrfl] m + ot _ m - 1, ) 

m>0 

These operators form a set of mutually commuting SU(2) algebras: 

(18) [S+, Sr] = 25ij§$, [Sf, Sf] = ±%5±. 
The operator S? can be related to the number operator 

(19) 5° = Nj hl h 

where Clj = j + \ is the maximum number of pairs that can occupy the level j and the 
number operator is 

(20) Nj = ^ (4 rrfij rn + a] _ m dj _ m ) . 

m>0 

Since < Nj < fij these SU(2) algebras are realized in the representation with the total 
angular momentum quantum number 5%. 

The most general Hamiltonian for nucleons interacting with a pairing force can be 
written as 

(21) # = E "A m*i m-\G\J2 CJJ> S+S-, , 

jm jj' 

where Cj are the single particle energy levels, \G\ is the pairing interaction strength with 
dimensions of energy, and Cjj> arc dimcnsionlcss parameters describing the distribution 
of this strength between different orbitals. When the latter are separable {cjji — c*Cj>) 
we get 

(22) H = J2 m *i m-\G\J2 ^i'S+Sj,. 

3™ jj' 
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There are a number of approximations one can make to simplify the Hamiltonians in Eqs. 
(21) or (22). If we assume that the NN interaction is determined by a single parameter 
(usually chosen to be the scattering length), all c/s are the same and we get 

(23) // (: Y, s ; s /- 

jm jj' 

This case was solved by Richardson [21]. 

If we assume that the energy levels are degenerate then the first term is a constant 
for a given fixed number of pairs. This case can be solved by using the quasispin algebra 
since H oc S + S~ . A list of the exactly solvable cases can then be given as follows: 

• Quasi-spin limit [20]: 

(24) H = -\G\J^SfST,. 

jj' 

• Richardson's limit [21]: 

(25) // (: Y, s ; s ,- 

jm jj' 

• An algebraic solution developed by Gaudin [22] which is sketched out in the next 
section (section 2"2). Gaudin's model is closely related to Richardson's limit. 

• The limit in which the energy levels are degenerate (when the first term becomes 
a constant for a given number of pairs): 

(26) H = -\G\Y / c* j c j ,S+S7 l . 

jj' 

Different aspects of this solution was worked out in Refs. [23, 24, 25]. 

• Most general separable case with two shells [26]. 

Clearly we can solve the pairing problem numerically in the quasispin basis [27]. Note 
that it is also possible to utilize the quasi-spin concept for mixed systems of bosons and 
fermions in a supersymmetric framework [28] . 

2'2. Gaudin Algebra. - To study the problem of interacting spins on a lattice Gaudin 
introduced a method [22] closely related to the Richardson's solution. Our presentation 
of Gaudin's method here follows that given in Ref. [29]. This method starts with three 
operators J ± (A), J°(A), parametrized by one parameter A, satisfying the commutation 
relations 



(27) 



A — /i 
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(28) [AX),JH,)}=± J±{X IZ J ^\ 
and 

(29) [J°(A),J°( M )] = [J ± (A),J ± ( M )]=0. 

The Lie algebra depicted above is referred to as the rational Gaudin algebra. A possible 
realization of this algebra is given by 

N OO N q± 

(30) jO iX) = J2^ and J±(A) = ]T — 

i=l i=l 

where e» are, in general, arbitrary parameters. In Eq. (30), instead of the quasi-spin 
algebra one can obviously use any mutually-commuting N SU(2) algebras (in fact, copies 
of any other algebra) . The operator 

(31) H(X) = J°(A)J°(A) + ij + (A)J-(A) + \ J- (A) J+(A) 

is not the Casimir operator of the Gaudin algebra, but such operators commute for 
different values of the parameter: 

(32) [H{X),H{fj,)]=0. 
Lowest weight vector, |0), is chosen to satisfy the conditions 

(33) J"(A)|0)=0, and J°(A)|0) = W(A)|0). 
Hence it is an eigenstate of the Hamiltonian given in Eq. (31): 

(34) tf(A)|0) = [W/(A) 2 -M/'(A)]|0). 

To find other eigenstates consider the state |£) = J + (£)|0) for an arbitrary complex 
number £. Since 

(35) [H(X), J+(0] = (J + (A) J°(0 - J+(0 J°(A)) , 

we conclude that if W(£) = 0, then J + (£)|0) is an eigenstate of H(X) with the eigenvalue 

(36) El{X )=[W(Xf-W'(X)]-2^. 

Gaudin showed that this approach can be generalized and a state of the form 

(37) ...,&>= J + (£i)J + fa) ■ ■ ■ J + (6)|o > 
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is an eigenvector of H(X) if the numbers £1,62, • • • ,£n € C satisfy the so-called Bethe 
Ansatz equations: 

n 

(38) W (U=J2r— T f ° r « =1 .2."-.n- 
Corresponding eigenvalue is 

(39) E n (X) = [W(Xf W'(X)} 2 £ ^'f^ . 

To make a connection to Richardson's solution we define so-called 7\L-operators as 

(40) \im (X-e k )H(X) =K k . 

In the realization of Eq. (30), these operators take the form 

(41) 7^ = -2]T^-^. 

Taking the limits \i — > e k first and A — > ej second in Eq. (32), one easily obtains 

(42) [H(X),n k ] = 0, [Kj,K k ] = 0. 
One can also prove the equalities 

(43) ^7^ = 0, 

i 

and 

(44) 5^6*7^ = -2 ^ Si -Sj. 

i i^j 

A careful examination of the Gaudin algebra in Eqs. (27), (28), and (29) indicates 
that not only the operators J(A), but also the operators J(A) + c satisfy this algebra for 
a constant c. In this case the conserved quantity is replaced by 

(45) H(X) = J(A) • J(A) => H(X) + 2c • J(A) + c 2 

which has the same eigenstates. One can define Richardson operators, R k , in an analo- 
gous way to the lZ k defined in Eq. (40): 

(46) lim(A-e fc )Gff(A)+2c-S)=i2 fc) 
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which implies 



(47) R k = -2c ■ S k - 2J2 
Eq. (42) is then replaced by 

(48) [H(X) + 2c • S, R k ] = [Rj , R k ] = 0, 
with the conditions 

(49) ^Ti? 4 = -2c.]Ts fe , 

i k 

and 

(50) ^e i fl i =-2^e i c.S i -2^S i -S j . 

i i ijtj 

Rewriting the Hamiltonian of Eq. (25) in the form 

(si) = £ *rf - w ( (e s *) • (e s ») - (E s?) 2 + (E s?) 

+ constant terms, 
and choosing the constant vector of Eq. (45) to be 
(52) c=(0,0,-l/2|G|) 
one immediately obtains 

< 53 » m 



= £e,R, + |G| 2 fe fi <) -l G E* + 

i \ i / i 



Since all R k and i/(A) + 2c • S mutually commute (cf. Eq. (48)), they have the same 
eigenvalues. Eq. (53) then tells us that they are also the eigenvalues of the Richardson 
Hamiltonian, Eq. (25). 

2'3. Exact Solution for Degenerate Spectra. - If the single particle spectrum is degen- 
erate, then it is possible to find the eigenvalues and eigenstates of the Hamiltonian in 
Eq. (22). If all the single particle energies are the same, the first term in Eq. (22) is a 
constant and can be ignored. Defining the operators 

(54) S+(0) = £c*S+ and S"(0) = 
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the Hamiltonian of Eq. (22) can be rewritten as 

(55) H = -\G\S + (0)S-(0) + constant. 

(The operators in Eq. (54) are defined with argument for reasons explained in the 
following discussion). In the 1970's Talmi showed that, under certain assumptions, a 
state of the form 

(56) £+(o)|o) = ]T c ;£+|o>, 

3 

with |0) being the particle vacuum, is an eigenstate of a class of Hamiltonians including 
the one above [30] . A direct calculation yields the eigenvalue equation 

(57) ff£+(0)|0) = ^-IGl^lc/j £+(0)|0>. 

However, there are other one-pair states besides the one in Eq. (57). For example for 
two levels j\ and j% , the orthogonal state 




is also an eigenstate with E=0. In Rcf. [23] it was shown that there is a systematic way 
to derive these states. To see their solution we define the operators 

(59) S + (a;)= V °? , 2 S+ and g»=^ °f Sr. 

^-r 1 1 — c,- rx J ' 1 — \Cjrx J 

3 3 

Note that if one substitutes x = in the operators of Eq. (59), one obtains the operators 
in Eq. (54). Further defining the operator 



one can prove the following commutation relations: 

(61) [S + (x),S-(0)} = [S+(0),S-(x)} = 2K°(x) 

(62) [K°(x),sHy)} = ± S±{x) ~ S±iy) 
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These commutators are very similar to, but not the same as, those of the Gaudin algebra 
described in Eqs. (27), (28), and (29). Using the commutators in Eqs. (61) and (62) one 
can easily show that 

(63) 5+(0)S+(^ JV) ).-^ + (4 A 2i)|0) 

is an eigenstate of the Hamiltonian in Eq. (55) if the following Bcthc ansatz equations 
are satisfied: 

-0-/2 1 N ~ 1 1 

(64) E 1/y2 E IW Jn) m = l,2,...N-l. 



The energy of the state in Eq. (63) is 
(65) E N = -\G 




The authors of Ref. [23] used a Laurent expansion of the operators given in Eq. (59) 
around x = followed by an analytic continuation argument to show the validity of their 
results in the entire complex plane except some singular points. The derivation sketched 
above instead utilizes the algebra depicted in Eqs. (61) and (62); it is significantly 
simpler. The state in Eq. (63) is an eigenstate if the shell is at most half full. Similarly 
the state 

(66) S+(4 N) )S^4 N) )...S+( X ^)\0) 

is an eigenstate with zero energy if the following Bethe ansatz equations are satisfied: 

(67) E 17| '^11*) = E Iw I Iw fOTeve T m = l,2,...,iV 

j 1 /\ C J\ x "i fe=l(fe^m) Xm X k 

Again this is an eigenstate if the shell is at most half full. 

To figure out what happens if the available states are more than half full we note 
that there are degeneracies in the spectra. Let us denote the state where all levels are 
completely filled by |0). It is easy to show that the state |0) and the particle vacuum, 
|0) are both eigenstates of the Hamiltonian of Eq. (55) with the same energy, E = 
— \G\ J2j Qj\ c j\ 2 - This suggests that if the shells are more than half full then one should 
start with a state of the form 

(68) S-^lS-^lJ-^llO). 
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Table I. - Particle-hole degeneracy - two different states with the same energy 



No. of Pairs 


State 


N 


s+( )s+(*r)---s + (4"-\)io> 




§-{zW)§-(4 tf) )...§-{zW 1 )\0) 



Indeed the state in Eq. (68) is an eigenstate of the Hamiltonian in Eq. (55) with the 
energy 

(69) E = -G (]T0,M 2 -£ 2 



(N) 

j fc=l *k 



if the following Bethe ansatz equations are satisfied [24, 25]: 



( 70 ) r A2_ AN) (N) + £ (JV)_ W 

j Vl^jl z "i fe=l(fe^m) Zm z fc 

In Eq. (70) N max + 1 — TV is the number of particle pairs. If the available particle states 
are more than half full, then there are no zero energy states. Note that the Bethe ansatz 
conditions in Eqs. (70) and (64) as well the energies in Eqs. (69) and (65) are identical. 
This particle-hole degeneracy is due to a hidden supersymmetry [25] and is illustrated in 
Table I. This supersymmetry is described in Section 3'2. 

2'4. Exact Solutions with two shells. - Consider the most general pairing Hamiltonian 
with only two shells: 

|§|=£ M° - £ <$«A *r + £ ^ 
j jj' j 

with £j = ej/\G\. It turns out that this problem can be solved using techniques illustrated 
in the previous sections. Eigcnstatcs of the Hamiltonian in Eq. (71) can be written using 
the new step operators [26] 

(72) = E o — %^r s t 

Y 2£j ~ ^ x 

as 

(73) J + (x 1 )J+(x 2 )---J+(x N )\0). 
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Introducing the definitions 



(74) /3 = 2 " J2 5 = 2 

we obtain the energy eigenvalue to be 

n=l 

if the parameters Xk satisfy the Bethe ansatz equations 

(76) " j|Cj|2 2 + E - • 

J — l y^pru ) 

2'5. Solutions of Bethe Ansatz equations. - Various solutions of the pairing problem 
discussed in the previous sections are typically considered as semi-analytical solutions 
since one still needs to find the solutions of the Bethe ansatz equations. This a task which, 
quite often, needs to be tackled numerically. However, in certain limits it is possible to 
find solutions of the Bethe ansatz equations analytically. The method outlined here was 
first presented in Ref. [31]. 

Consider the Bethe ansatz equations for the degenerate single particle levels and zero 
energy eigenstate given in Eq. (67). Introducing new variables, r]i, 

(77) 4 N) ~ 1 

I L J2 I 

Eq. (67) can be rewritten as 



' 78 ' ~Jn) Tn) TnY + " (jvJ — °- 

fe=W0»7i ~% Vi 1 — »7i 

It can be easily shown that the polynomial admitting the solutions of Eq. (78) as zeros 

N 

(79) Pn(z) = ]J(z-4 N) ) 

i=l 

satisfies the hypergeometric equation [32] 

(80) z(l - z)p% + [-% + (%%) z}p' N + N(N- n n - n J2 -l)p N = 0. 

Consequently the problem of finding the solutions of the Bethe ansatz equation (67) 
reduces to calculating the roots of hypergeometric functions. For analytical expressions 
of the energy eigenvalues obtained in this manner the reader is referred to Ref [24] . 
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3. Supersymmetric Quantum Mechanics in Nuclear Physics 

Consider two Hamiltonians 

(81) Hi = G f G, H 2 = GG\ 

where G is an arbitrary operator. The eigenvalues ol these two Hamiltonians 

(82) GtG|l,n)=£«|l,n>, 

Gtf\2,n) =E^\2,n) 

are the same: 

(83) EW = = E n 
and the eigenvectors are related: 

(84) |2,n) = G[GtG]" 1/2 |l,n). 

This works for all cases except when G|l,n) = 0, which should be the ground state 
energy of the positive-definite Hamiltonian Hi . 

To see why this is called supersymmetry we define 



Then 

(86) 
with 



Gt J ' * 10 



(87) [H,Q]=0=[H,Ql]. 

Clearly the operators H, Q, and close under the commutation and anticommutation 
relations of Eq. (86) and (87), forming a very simple superalgebra. The two Hamiltonians 
depicted in Eq. (81) are said to form a system of supersymmetric quantum mechanics 
[33]. Many aspects of the supersymmetric quantum mechanics has been investigated in 
detail [34, 35, 36]. In the following sections two applications of supersymmetric quantum 
mechanics to nuclear physics are summarized. 
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3*1. Application of Super symmetric Quantum Mechanics to Pseudo- Orbital Angular 
Momentum and Pseudospin. - The nuclear shell model is a mean-field theory where the 
single particle levels can be taken as those of a three-dimensional harmonic oscillator 
(labeled with SU(3) quantum numbers) for the lowest (A < 20) levels. For heavier 
nuclei with more than 20 protons or neutrons, different parity orbitals mix. The Nilsson 
Hamiltonian of the spherical shell model is [37] 

(88) H = ub\bi - 2/cL.S - fc^L 2 , 

where the second term mixes opposite parity orbitals and the last term mocks up the 
deeper potential felt by the nucleons as L increases. 

Fits to data indicate that [i rts 0.5, hence there exists degeneracies in the single 
particle spectra. In the 50-82 shell (the SU(3) label or the principal harmonic oscillator 
quantum number of which is N = 4) the S]/2 and d 3 / 2 orbitals and further d 5 / 2 and 
g7/2 orbitals are almost degenerate. It is possible to give a phenomenological account of 
this degeneracy by introducing a second SU(3) algebra called the pseudo-SU(3) [38, 39]. 
Assuming that those orbitals belong to the N — 3 (with I = 1,3) representation of the 
latter SU(3) algebra, the quantum numbers of the SO (3) algebra included in this new 
SU(3) are called pseudo-orbital-angular momentum [l = 1,3 in this case). We can also 
introduce a "pseudo-spin" (s = |). One can easily show that j = 1/2 and 3/2 orbitals 
(and also j = 5/2 and 7/2 orbitals) are degenerate if pseudo-orbital angular momentum 
and pseudo-spin coupling vanishes. This degeneracy follows form the supersymmetric 
quantum mechanical nature of the problem. 

It can be shown that two Hamiltonians written in the SU(3) and the pseudo-SU(3) 
bases are supersymmetric partners of each other [40] . The operator that transforms these 
two bases into one another is [41] 

(89) U = G [&G] ~ V2 = V2F_ {K + [F+, F_])" 1/2 

= (Wi) (bibi-o-iL z y 1/2 

yielding the supersymmetry transformation between the pseudo-SU(3) Hamiltonian H' 
and the SU(3) Hamiltonian H: 

(90) H' = U HU^ = b\bi — 2k (2/j, — 1) L • S — k/^L 2 + [1 - 2fc( M - 1)] . 

The transformation depicted in Eq. (90) can also be expressed in terms of the generators 
of the orthosymplectic superalgebra Osp(l/2) [40]. 

3'2. Supersymmetric Quantum Mechanics and Pairing in Nuclei. - Let us consider 
the separable pairing Hamiltonian with degenerate single-particle spectra given in Eq. 
(55): 



(91) 



Hsc~-\G\S+(0)S-(0), 
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and introduce the operator 




This operator transforms the empty shell, |0), to the fully occupied shell, |0): 

(93) f |0> = |0). 

To establish the connection to the supersymmetric quantum mechanics we define the 
operators 

(94) B-=f^S-(0), B+ = S+(0)f. 

Supersymmetric quantum mechanics tells us that the partner Hamiltonians Hi = B + B~ 
and H2 = B~B + have identical spectra except for the ground state of Hi. It can easily 
be shown that in this case two Hamiltonians Hi and H2 are actually identical and equal 
to the pairing Hamiltonian (Eq. (91)). Hence the role of the supersymmetry is to connect 
the states |l,n) and \2,n). These are the "particle" and "hole" states. Hence if the shell 
is less than half-full ( "particles" ) there is a zero-energy state (note that the Hamiltonian 
in Eq. (55) negative definite, hence this is the highest energy state). If the shell is more 
than half- full ("holes") this state disappears. Otherwise the spectra for particle and hole 
states are the same as the rules of the supersymmetric quantum mechanics implies [25]. 

4. Dynamical Supersymmetries in Nuclear Physics 

Dynamical supersymmetries in nuclear physics start with the algebraic model of 
nuclear collectivity called Interacting Boson Model [6, 7, 8]. In this model, low-lying 
quadrupole collective states of even-even nuclei are generated as states of a system of 
bosons occupying two levels, one with angular momentum zero (s-boson) and one with 
angular momentum two (d-boson). It is discussed elsewhere in great detail in these pro- 
ceedings [42, 43]. There are three exactly solvable limits of the simplest form of the 
Interacting Boson Model where neutron and proton bosons are not distinguished: 

• Vibrational Limit: SU(6) D 5t/(5) D 50(5) D 50(3). 

• Rotational Limit: SU(6) D SU(3) D 50(3). 

• Gamma-Unstable Limit: SU(6) D 50(6) D 50(5) D 50(3). 

It is possible to extend the Interacting Boson Model to describe odd-even and odd-odd 
nuclei [44] . Dynamical supersymmetries arise out of certain exactly solvable limits of this 
Interacting Boson- Fcrmion Model [45, 46]. 

In an odd-even nucleus, in addition to the correlated nucleon pairs (s and d bosons) 
we need the degrees of freedom of the unpaired fcrmions. If those unpaired fermions are 
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in the ji, J2, J3, • • • orbitals then the fermionic sector of the theory is represented by the 
fcrmionic algebra <S'?7p'(^ i (2j i + l)) and the resulting S(7(6)sX SUF(%2i(2ji + l)) algebra 
is embedded in the superalgebra 51/(6/ ^2A2ji + 1))- In the first example of dynamical 
supersymmetry worked out the unpaired fermion was in a j = 3/2 (d 3 / 2 ) orbital coupled 
to the nuclei described by the gamma-unstable (SO (6)) limit of the Interacting Boson 
Model. The resulting S?7(6/4) supersymmetry was used to describe many properties 
of nuclei in the Os-Pt region [47]. Immediately afterwards this supersymmetry was 
extended to the SU (6/ 12) superalgebra including fcrmions in S1/2, d.3/2, and d.5/2 orbitals 
[48]. Theoretical implications of nuclear supersymmetries were extensively investigated 
by many authors [49, 50, 51, 52, 53, 54]. The existence of dynamical supersymmetries in 
nuclei is experimentally established [55, 56, 57, 58, 59, 60]. 

There exists other symmetries of nuclei describing phase transitions between different 
dynamical symmetry limits. For these so-called critical point symmetries zeros of wave- 
functions in confining potentials of the geometric model of nuclei give rise to point groups 
symmetries [2] . This scheme is generally applicable to the spectra of systems undergoing 
a second-order phase transition between the dynamical symmetry limits SU(n — 1) and 
SO(n). The resulting symmetries are either named after the discrete subgroups of the 
Euclidean group E(5), or X(5). A review of the experimental searches for critical point 
symmetries at nuclei in the A=130 and A=150 regions is given in Ref. [61]. It is also pos- 
sible to couple Bohr Hamiltonian of the geometric model with a five-dimensional square 
well potential to a fermion using the five-dimensional generalization of the spin-orbit 
interaction. In doing so the E(5) symmetry of the even-even nuclei goes into the E(5/4) 
supersymmetry for odd-even nuclei near the critical point [62] . Initial experimental tests 
of this latter supersymmetry are encouraging [63] 

5. Application of Symmetry Techniques to Subbarrier Fusion 

In the study of nuclear reactions one needs to describe translational motion cou- 
pled with internal degrees of freedom representing the structure of colliding nuclei. The 
Hamiltonian for such a multidimensional quantum problem is taken to be 

h 2 

(95) H = -— V 2 + V(r) + H (q) + H int (r, q), 

where r is the relative coordinate of the target-projectile pair and q represents any internal 
degrees of freedom of the system. To study the effect of the structure of the target nuclei 
on fusion cross sections near and below the Coulomb barrier one can take V(r) to be 
the potential barrier and the third term in Eq. (95) represents the internal structure of 
the target nuclei. All the dynamical information about the system can be obtained by 
solving the evolution equation 



(96) 
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with the initial condition U(ti) = 1. If part of the Hamiltonian in Eq. (95), say 
H {q) + £/i n t(r, q), is an element of a particular Lie algebra, then symmetry methods 
can significantly simplify the solution of the problem. In such a case, if one considers a 
single trajectory r(t), the q-dependent part of the evolution operator becomes an clement 
of the Lie group associated with the Lie Algebra mentioned above. (It is possible to make 
this statement rigorous in the context of path- integral formalism [64] ) . 

It is by now well-established that heavy-ion fusion cross-sections below the Coulomb 
barrier are several orders of magnitude larger than one would expect from a one dimen- 
sional barrier penetration picture, an enhancement which is attributed to the coupling of 
the translational motion to additional degrees of freedom such as nuclear and Coulomb 
excitation, nuclcon transfer, or neck formation [65, 66]. The multidimensional barrier 
penetration problem inherent in subbarrier fusion can be addressed in the coupled- 
channels formalism and state of the art coupled-channel codes are currently available 
(see e.g. [67]). Although several puzzles remain (such as the large values of the surface 
diffuseness parameter in the nuclear potential required to fit the data [68, 69] ; very steep 
fall-off of the fusion data for some systems at extreme subbarrier energies [70] ; or inade- 
quacy of the standard nuclear potentials to simultaneously reproduce fusion and elastic 
scattering measurements [71]) there is overall good agreement between coupled-channels 
calculations and the experimental data. 

An alternative approach is to formulate the problem algebraically by using a model 
amenable to such an approach for describing the nuclear structure, such as the Interacting 
Boson Model [72] . Not only this approach fits the data well, it can also be used for tran- 
sitional nuclei, the treatment of which could be more complicated for coupled-channels 
calculations [72, 73, 74, 75]. 

6. Application of Algebraic Techniques in Nuclear Astrophysics 

Nuclear astrophysics has been very successful exploring the origin of elements. As our 
understanding of the heavens evolved, it was realized that nuclear data was needed for 
astrophysics calculations, such as nucleosynthesis, stellar evolution, the Big Bang cosmol- 
ogy, x-ray bursts, and supernova dynamics. Nuclear astrophysics initially started with 
reaction rate measurements, but with the recent rapid growth of the observational data, 
expanding computational capabilities, and availability of exotic nuclear beams, interest 
in all aspects of nuclear physics relevant to astrophysical phenomena has significantly 
increased [76]. In the rest of this section an application of algebraic techniques to a 
nuclear astrophysics problem is presented. 

Light nuclei are formed during the big-bang nucleosynthesis era and nuclei up to the 
iron, nickel and cobalt group are formed during the stellar evolution. A good fraction 
of the nuclei heavier than iron were formed in the rapid neutron capture (r-process) 
nucleosynthesis. The astrophysical location of the r-process is expected to be where 
explosive phenomena are present since a large number of interactions are required to 
take place at this location during a rather short time interval. Core-collapse supernovae 
which occur following the stages of nuclear burning during stellar evolution after the 
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formation of an iron core is such a site. Neutrino interactions play a very important role 
in the evolution of core-collapse supernovae [77]. Almost all (99%) of the gravitational 
binding energy (10 53 ergs) of the progenitor star is released in the neutrino cooling of 
the neutron star formed after the collapse. It was suggested that neutrino-neutrino 
interactions could play a potentially very significant role in core-collapse supernovae [78, 
79, 80, 81, 82]. A key quantity for determining the r-process yields is the neutron to seed- 
nucleus ratio (or equivalently neutron-to-proton ratio). Interactions of the neutrinos and 
antineutrinos streaming out of the core both with nucleons and seed nuclei determine the 
neutron-to-proton ratio. Before these neutrinos reach the r-process region they undergo 
matter-enhanced neutrino oscillations as well as coherently scatter over other neutrinos. 
Many-body behavior of this neutrino gas is not completely understood, but may have a 
significant impact on r-process nucleosynthesis. 

It was shown that neutrinos moving in matter gain an effective potential due to 
forward scattering from the background particles such as the electrons [83, 84, 85]. One 
can write the Hamiltonian describing neutrino transport in dense matter for two neutrino 
flavors as 



where N e is the electron density of the medium (assumed to be charge neutral and 
unpolarized), 9 is the mixing angle between neutrino flavors, Sm 2 — m\ — m\ is the 
difference between squares of the neutrino masses, and Gf is the Fermi coupling strength 
of the weak interactions. In the above equation the operators J ± , J° has been written 
in terms of the neutrino creation and annihilation operators: 



Here a\(p) and a x (p) are the creation and annihilation operators for the electron neutrino 
with momentum p and either muon or tau neutrino with momentum p, respectively. The 
operators in Eq. (98) form as many mutually commuting SU(2) algebras as the number 
of allowed values of neutrino momenta: 



If, instead of two neutrino flavors, one considers all three active neutrino flavors then 
one should use copies of SU(3) algebras. The contribution of neutrino-neutrino forward 
scattering terms to the neutrino Hamiltonian is given by 



(97) 




(98) 



J+(p) = al(p)a e (p), J-{p) = al(p)a x (p) 
Mp) = g (4(p)°*(p) - al{p)a e (p)) . 



(99) 



[J+(p), J_(g)] = 2S 3 (p - q)J (p), [J (p), J ± (p)] = ±<5 3 



(p-q)J±(p) 



(100) 
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where d pq is the angle between neutrino momenta p and q. It is easy to include neutrino- 
antineutrino and antineutrino-antineutrino scattering terms in this Hamiltonian [81], but 
we ignore them here to keep the presentation simple. 

Since a large number of neutrinos (10 58 ) are emitted during a typical core-collapse it 
is very difficult to evaluate neutrino evolution exactly using the two-body Hamiltonian 
of Eq. (100). Instead one uses a mean field approximation where the product of two 
commuting arbitrary operators 0\ and (D 2 can be approximated as 

(101) <9A ~ d.imo + (mods - miomio, 

provided that the condition 

(102) (^\did 2 \o = (momo 

is satisfied. In Eqs. (101) and (102) One can then replace the exact Hamiltonian of Eq. 
(100) with the approximate expression 

H vv ~ 2^^ J dpdqR pq (j (p)(Mq)) + l -J + {p)(J-{q)) + ±J_(p) {./+(«)>) 

where R pq = (1 — cos, , d pq ) and the averages are calculated over the entire ensemble of 
neutrinos. Systematic corrections to this expression are explored in [81]. Calculations 
using this approximation do not yield conditions (i.e. large enough neutron to seed 
nucleus ratio) favorable to r-process nucleosynthesis [78]. There is encouraging progress 
in numerical calculations using the exact Hamiltonian of Eq. (100) [79, 80]. However an 
algebraic solution to this problem is currently lacking. 

7. Conclusions 

In many-fermion physics, mean field approaches can describe many properties; but 
are inadequate to describe the whole picture; pairing correlations play a crucial role. In 
fact, pairing correlations are not only necessary to understand the structure of rare-earths 
and actinides, but, since they arc essential for the description of the neutrino gas in a 
core-collapse supernova where many nuclei are produced, also necessary to understand 
the existence of such nuclei in the first place. Models exploiting symmetry properties 
and pairing correlations have been very successful. These also gave rise to dynamical 
supersymmetries. We showed that using algebraic techniques it is possible to solve the 
s-wave pairing problem almost exactly (i.e., reducing it to Bethe ansatz equations) at 
least for a number of simplified cases. 

* * * 
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